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This is a doctorate level lecture on the physics of accretion discs driving 
magnetically self-confined jets, usually referred to in the literature as disc 
winds. 

I will first review the governing magnetohydrodynamic equations and then 
discuss their physical content. At that level, necessary conditions to drive jets 
from keplerian accretion discs can already be derived. These conditions are 
validated with self-similar calculations of accretion-ejection structures. 

In a second part, I will critically discuss the biases introduced when using 
self-similarity as well as some other questions such as: Are these systems really 
unstable? Can a standard accretion disc provide the conditions to launch jets 
in its innermost parts? What is the difference between X-winds and disc- 
winds? 

Finally, the magnetic interaction between a protostar and its circumstellar 
disc will be discussed with a focus on stellar spin down. 



1 The accretion-ejection paradigm 

This lecture is designed to be read with an accompanying file (pdf or ppt) 
where more illustrations and figures can be found. It can be retrieved at the 
URL: http://www-laog.obs. ujf-grenoble.fr/^ferreAra/JETSET/ 'school.html I 
also recommend the reviews of Konigl & Pudritz 2000; Ferreira 2002. In 
iFerreira et alJ l)2006a(l . a review of all MHD models for Young Stellar Objects 
has been made with a comparison of the corresponding jet kinetic observa- 
tional properties. 



1.1 A "universal" picture 

Actively accreting "classical" T Tauri stars (TTS) often display supersonic 
collimated jets on scales of a few 10-100 AU in low excitation optical forbid- 
den lines. Molecular outflows observed in younger Class and I sources may 
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be powered by an inner unobserved "optical jet" (see Cabrit's contribution, 
this volume). These jet signatures are co rrelated with t he infrared excess and 
accre tion rate of the circumstellar disc jCabrit et all Il990t lHartigan et all 
19951) . It is therefore widely believed that the accretion process is essential to 
the observed jets, although the precise physical connection remains a matter 
of debate: do the jets emanate from the star, the circumstellar disc or the 
magnetospheric star-disc interaction? 

One argumen t in favor of accretion-powered disc winds is its "universal- 
ity" l)Liviol 1 1997^1 . Indeed, self-collimated jet production from accretion discs is 
also invoked to explain an accretion-ejection correlation observed in compact 
objects (i.e. some a ctive galactic nuclei, quasars and X-ray binaries, see eg. 
iMerloni et al.ll2003l and references therein) . The underlying idea is quite sim- 
ple: accretion discs around a central object can, under certain circumstances 
and whatever the nature of this object, drive jets through the action of large 
scale magnetic fields. These fields would tap the mechanical energy released 
by the mass accreting in th e disc and transfer it to the fraction that is ejected 
i Blandford fc PavneL Il982l) . The smaller the fraction and the larger the final 
jet velocity. One thing that must be understood is how the presence of such 
jets modifies the nature of the underlying accretion flow. Many papers in the 
literature actually assume (implicitly or not) th at the accretion disc resem - 
bles a standard ac cretion disc as envisioned bv IShakura fc Sunvaevl l)l973|) : 
iFrank et al.l ((2002) . This is wrong as will be shown later. 

A Magnetized Accretion-Ejection Structure (hereafter MAES) is an accre- 
tion disc where accretion and ejection are interdependent processes. As such, 
it is composed of an accretion disc (called hereafter JED for Jet Emitting 
Disc) thread by a large scale magnetic field of bipolar topology and giving 
rise to the two bipolar jets. The goal of the study of a MAES is to obtain 

(1) - the conditions allowing for a steady state accretion-ejection process; 

(2) - the ejection to accretion rates ratio as function of the disc physical con- 
ditions; 

(3) - the jet properties (kinematics, power, shape) as function of the disc prop- 
erties. 

1.2 From magnetostatics to magnetohydrodynamics 

Magnetohydrodynamics (MHD) is the theoretical framework required to de- 
scribe the interaction between an ionized gas and magnetic fields. But mag- 
netostatics is very helpful to understand basic mechanisms. 

A zeroth order description of a MAES is that of a rotating conducting 
disc thread by a magnetic field aligned with the rotation axis (much alike 
a Barlow's wheel). According to Faraday's induction law, an electromotive 
force (emf) across the disc, e = J(uA B) ■ dr — J QrB z dr, creates an electric 
potential difference between the disc center and its border (Fig. QJ. If some 
conducting wire connects the border to the center, closing thereby an electric 
circuit, then a radial electric current is induced. Because of this current I, the 
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disc becomes prone to a Laplace force, F = J IB z dr, which will slow down 
the disc (Lenz's law). One could also say that the field "resists" to the shear 
provoked by the rotation (the current / induces a toroidal component B^). But 
such a " mechanical" view of the magnetic field disregards its electromagnetic 
nature and one may tend to forget that electric currents must be maintained 
and able to flow... 




V u 



Fig. 1. Left: A rotating disc embedded in a magnetic field induces 
a current leading to a magnetic braking (Barlow's wheel: see e.g. 
http://www.sparkmuseum.com/MOTORS.HTM for many illustrations). Right: 
A MAES can be seen as two independent electric circuits, each corresponding to a 
jet. Asymmetric jets can thus be easily achieved, even with a symmetric poloidal 
field. 



In astrophysics, the disc is made of gas that, provided it can cross the field 
lines, will accrete towards the central object as it looses angular momentum. 
This angular momentum is linked to the electric current flowing in the jets: 
the jet kinetic power is fed by the flux of magnetic energy provided at the disc 
surface. Note that while the streamlines of the ejected material go to infinity 
those of the current density must be closed and return to the disc where the 

emf is. 

T his is actually the reaso n why jet collimation is a subtle issue ijHevvaerts fc Normal 

Il989t l20o3 lOkamotol l2003|b Make a cut at a distance z of a jet and compute 
the total current flowing inside it, namely / = J dr2irr J z . If this current is 
non zero and (for instance) negative, then one might say that the Laplace 
force will be directed towards the jet axis (Ampere's theorem tells that is 
negative in that case). This is the basic idea of the "magnetic hoop stress" 
that provides a self-confinement to jets. However, the local magnetic force is 
actually J A B and depends on the radial distribution of J z (r) ! This depends 
on the lateral boundary conditions (jet axis and outer edge) but also on what 
happened upstream (or in the past, if we follow a lagrangean particle): since 
jet acceleration is a conversion of electric into kinetic power, then jet colli- 
mation depends as much on jet acceleration. One cannot therefore solve the 



4 



Jonathan Ferreira 



jet problem assuming for instance the shape of the field lines: the full MHD 
equations must be solved. 

The current flowing inside jets is precisely the current that allows for 
accretion. The accretion-ejection phenomenon has therefore to be viewed as 
a global electric current system. 

1.3 Basic assumptions 

Modeling a MAES requires several assumptions: 

(1) Presence of a large scale vertical magnetic field in the disc. Its 
origin and amplitude remain an open question. For the purpose of illustration, 
we will assume a positive vertical component B z anchored in the disc (bipolar 
topology). 

(2) Single-fluid MHD: matter is assumed ionized enough so that all 
species (ions, neutrals and electrons) are well coupled and can be treated as a 
single fluid. Such an assumption shou ld always be verified a posteriori for any 
model but it is seldom made (see e.g. lGarcia et al.|l200ll for how to do it). 

(3) Axisymmetry: using cylindrical coordinates (r, 0, z) all quantities 
are assumed to be independent on <j), the jet axis being the vertical axis. Then, 
Ecj, — and all quantities can be decomposed into poloidal (the (r, z) plane) 
and toroidal components, eg. u = u p + flre^ and B — B p + B$e$. A bipolar 
magnetic configuration can then be described with B p = ^Va A e^, where 
the magnetic flux function a{r, z) is an even function of z and with an odd 
toroidal field B^(r, -z) = -B$(r, z) (Fig. 01. 

(4) Non-relativistic MHD, since observed motions are non-relativistic 
(this criterion is enough as long as MHD ordering applies). 

(5) Steady-state: all astrophysical jets display proper motions and/or 
emission nodules, showing that they are either prone to some instabilities or 
that ejection is an intermittent process. However, the time scales involved in 
all objects (from 1 to 10 2 yrs) are always larger than the orbit time scales 
in the innermost regions of the underlying accretion disc (close to the star). 
Therefore, a steady state approach is appropriate as a first step, while numer- 
ical simulations will be required to investigate time-dependent flows. 

1.4 Governing MHD equations 

According to the aforementioned assumptions, we use the following set of 
MHD equations (in MKSA units): 
Mass conservation 



V ■ pu = 



(1) 



Momentum conservation 



pu ■ Vm 



VP - pV$ G + J A B + V • T 



(2) 
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Fig. 2. Axisymmetric jets are made of magnetic surfaces of constant magnetic flux 
nested around each other and anchored in the disc. Each surface behaves like a 
funnel whose shape depends on the transfield equilibrium. Solving the jet equations 
requires to specify several quantities (see text). 

where <!?g is the central star gravitational potential and J = V A B/fi Q is the 
electric current density. The last term (with the stress tensor T) is actually due 
to a sustained turbulence inside the disc (it vanishes outside) which allows to 
transport angular momentum radially in the outward direction (see Terquem's 
contribution). It is presumably due to the presence of small scale magnetic 
fields but is usually grossly modeled by an anomalous viscosity v v — a v C s h, 
where a v is a free pa rameter, C K the disc sound speed and h ( r) the local disc 
vertical scale height ((Shakura fc Sunvaevi Il973 iFrank et all 120021) . 
Ohm's law and toroidal field induction 1 



where r\ m = [i, Vm and r( m = [i Q v' m are anomalous magnetic resistivities. The 
origin of these resistivities is the same as for viscosity, namely turbulence and 
they also vanish outside the disc. One expects turbulent media to display 
anomalous transport effects of heat, momentum but also magnetic flux. Note 
however that rotation in a Keplerian accretion disc introduces a strong dy- 
namical constraint. Indeed, the shear induced by rotation will unavoidably 
lead to huge toroidal magnetic fields until reconnection takes place (triggered 
by e.g. the tearing mode instability). As a consequence the amount of mag- 
netic dissipation in the toroidal direction might be much larger than in the 
poloidal direction. This has lead to the introduction of two anomalous coeffi- 
cients, ?2HLli^mJ_2ii^J2mJ£pvLi!£i a t e d respectively to the poloidal and toroidal 
fields llFerreira & Pclletier. 1993). 




Perfect gas law 



1 See Pelletier's contribution in this volume. Remember that E^, = while some 
algebra is required in order to derive Eq.^Jfrom the induction equation. 



VmJ<p = u p A B p 



(3) 



V • (%VrB ) = V • -{B^Up - Bp fir) 



(4) 
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P = p 



k B 



T 



(5) 



fim p 



where m p is the proton mass and p, a generalized "mean molecular weight" 
(in a fully ionized plasma p = 1/2). This expression assumes that all flu- 
ids (electrons, neutrals and ions) have the same temperature T. This is ful- 
filled only if the thermalization time scale (usually done through collisions) 
is short enough. Su ch an assumption should always be verified a posteriori 



where Q = Q + — Q~ is the sum of all heating Q + and cooling Q~ terms 
(including thermal conduction) and 7 the adiabatic index. There are many 
unsolved issues related to this exact equation for a single fluid. 

(1) Inside the disc, turbulence leads to an energy dissipation Q + = rj m J 2 + 
rj' m Jp+pu v \r\/ 'f2\ 2 , respectively Joule and "viscous" heating, but also to a cool- 
ing due to an energy transport by anomalous thermal conductivity. Moreover, 
the disc being optically thick, the radiation transport critically depends on the 
local opacity regime, which varies both with radius and height. Moreover, the 
disc surface is also the opticall y thick-thin transition, which is always an issue 
( see iFerreira fc Peiletierlll995l for a discussion). Besides, the energy equation 
in a standard accretion disc is usually writ ten Q + = Q~, the other terms 
being of the order (h/r) 2 ijFrank et all 120021. But these terms are important 
in the jet and cannot be neglected. 

(2) In the jet itself, although radiation may not be the dominant cooling term, 
it must be taken into account if one desires to compute e.g. the jet (forbid- 
den or permitted) emission lines or even radio continuum. A realistic and 
self-consistent treatment of the energy equation is therefore still out of range 
(even if one decouples the disc and its jets) and some stratagems must be 
used. 

The simplest way to deal with the energy equation in a MAES (valid in 
both the disc and its jets) is to use a polytropic equation of state P = Kp r , 
where the polytropic index r can be set to vary between 1 (isothermal case) 
and 7 (adiabatic case) . Note that K has to vary radially but remains constant 
along each field line: the jet entropy is thus fixed by the conditions prevailing 
at the disc surface. 

A more sophisticated approach can be done by prescribing the function 
Q along the field lines (this is equivalent to prescribing a variation of the 
polytropic index r). This will be discussed further in Section 3. 




Energy equation 



T 

Pu • V In — = (7 - 1) (Q + Pu • Vlnp) 



(G) 
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2 Physics of Jet Emitting Discs 

In this section, I will briefly discuss all the relevant physical effects that have 
to be covered in order to consistently describe Jet Emitting Discs or JEDs. 




Fig. 3. Sketch of the Jet Emitting Disc (JED) established between n and r e . The 
surface of the jet is determined by the magnetic surface anchored on r e . While the in- 
ner radius n is probably defined by some equilibrium with the stellar magnetosphere 
(see Section 4.6), the outer radius r e is free (it depends mostly on the magnetic flux 
available in the disc) . 



2.1 Mass conservation 

The disc accretion rate is defined as M a = —2 J Q 2irrpu r dz. In a standard 
accretion disc (hereafter SAD) M a is a constant both in time and radius. 
On the contrary, a JED displays mass loss at its surfaces so that M a must 
vary with the radius. This mass loss is parametrized by M (r) oc where the 
ejection index £ > is a measure of the disc ejection efficiency: the larger £ the 
larger the mass loss. The global mass conservation in a JED is M a (r e ) — 2Mj = 
M a (ri) where r e and are respectively the outer and inner radii of the JED 
and Mj is the mass flux from one side of the disc. The ejection to accretion 
rates ratio is 2Mj/M a (r e ) ~ £ln(rj/r^) and depends on both the ejection 
index £ and the radial extent of the JED (it will be shown later that £ is 
smaller than unity). The goal is of course to compute £ as a function of the 
disc physical conditions. 

2.2 Poloidal field diffusion 

Let us assume a smooth flux function a (see Sect. 1.3) so that a(r, z) ~ 
o (r)(l — z 2 /2l 2 ) where l(r) is the magnetic flux vertical scale height. Then, the 
bending of the poloidal field lines is measured at the equatorial plane by the 
magnetic Reynolds number lZ m = —ru r /is m = r 2 /I 2 . Such a bending is due to 
the interplay between advection by the accreting material and the turbulent 
magnetic diffusivity v m . It has been prescribed with v m = a m V/Ji, where 
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Va = B z / \fHoP is the Alfven speed at the disc midplane I Ferreira fc Pelletieil 
llflfl.^ . 

Now, magneto-centrifugal accel eration requires field lines bent enough at 
the disc surface, namely > B z l|Blandford fc Pavnelll982l quantities eval- 
uated at z = h are denoted with a superscript "+"). Since B+ /B z ~ H m h/r 
this implies lZ m > r/h. 

2.3 Angular momentum conservation 

The disc angular momentum can be transported by two means: (a) radially 
through a "viscous" turbulent torque which is probably triggered and sus- 
tained by an MHD insta bility such as the magneto-rotational instability (see 
Terquem's contribution. iBalbusI 120031 and references therein); (b) vertically 
by the jets. The viscous torque writes F V i SC ^ ~ —a v P/r where P is the total 
(gas+radiation) pressure and a v the so-called Shakura-Sunyaev parameter. 
The torque due to the jets writes F ma g,4> = JzB r — J r B z and its vertical be- 
havior strongly depends on the radial current density J r . At the disc midplane 
Fmag.cf) = ~JrB z ~ B~f B z j \x h and the disc angular momentum conservation 
reads 

1 + A ^^ = n e = n m ('^-) (7) 

where v v is the turbulent "viscosity" and A — F m ag,<j>/ Fvisc,<p- In a turbulent 
medium, one usually assumes that all anomalous transport coefficients are of 
the same magnitude so that v m ~ v v . In that situation, one gets the following 
consequences: 

- In a SAD, there is no jets and A = 0. Then lZ m ~ TZ e ~ 1 and, indeed, fiel d 
lines are too straight for a magneto-centrifugal driving l|Lubow et allll994a|k 

- In a JED, jets require lZ m ~ lZ e > r/h and thus A ~ r/h 3> 1: all the 
angular momentum must then be carried away by the jets, which results in 
an accretion velocity much larger than in a SAD. The "vi scous" torque is 
totall y negligible (in contrast to what is often assumed, e.g. lOeilvie fc Livid 
1998). 

This very important constraint (A > r/h ) can only be achieved if 
—B^B z /[i ~ P, that is with equipartition fields l|Ferreira fc Pelletierl Il995^ . 
Let us introduce here two important parameters: the disc magnetization 
[i = B^/fj, P and the magnetic shear q ~ —B^/B z . If q^i is not close to 
unity then no magncto-ccntrifugally driven jets can be launched from accre- 
tion discs. 



2.4 Toroidal field induction 

Magnetic driving of jets requires that the magnetic field starts to accelerate 
material at the disc surface. Hence, a JED must provide a transition from 
F m ag,<t> < at z = to F mag<< f, > at z = h and beyond. The only way to 
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achieve this is by allowing J r to decrease on a disc scale height. The vertical 
profile of J r is pr ovided by the induction e quation Q which, in a thin accretion 
disc, writes f see IFerreira fc Pelletierlll995|) 

r/' rn J r ~ i]' J + r [ dzB p -Vn- B 4> u z (8) 
Jo 

where rj' J = r)' m J r (z = 0). With no differential rotation, J r would remain 
constant and so would F mag ^ (< 0). In order to make J r decrease on a disc 
scale height, the disc differential rotation term must balance the term rf a J , 
due to the Faraday's induction law (the Barlow's wheel current). This can 
be done as long as —B^/B+ ~ l/a m (with u m ~ v' m ). Using the fact that 
.B+ > B z , one gets a magnetic shear q > l/a m . 



2.5 Disc vertical equilibrium 

It is worthwhile to consider the following general equality 

(J A B p ) ■ B^ = — (J p A B^) ■ B p (9) 

When the magnetic torque (lhs) is negative, so must be the projection of 
the Lorentz force on the poloidal field (rhs). Thus, deep within the disc, the 
poloidal Lorentz force is directed outwardly and towards the disc midplane. 
A quasi MHS equilibrium is therefore established with the balance between 
the total (gas+radiation) pressure gradient on one side and the magnetic 
compression due to the radial and toroidal field components and the gravity 
on the other side. Now, as one goes up in z and the magnetic torque changes its 
sign, the disc material starts to be azimuthaly accelerated. Correspondingly, 
the projection of the Lorentz force becomes also positive and helps to lift 
material out of the disc. 

This can be done in two ways (Fig. |U): (a) with a negative vertical com- 
ponent of the Lorentz force but a large radial component; (b) with a positive 
vertical component and a smaller negative radial component. Case (a) cor- 
responds to a small mass flux (£ < 1/2) where disc material must be lifted 
against the magnetic compression by the sole effect of the (gas+radiation) 
pressure gradient. Case (b) leads to a large mass flux (£ > 1/2) because of 
the magnetic pull due to the toroidal field pressure. 

In fact, it can be shown analytically that only solutions with £ < 1/2 can be 
stationary: solutions wit h large mass flu xes do not have enough power to allow 
for super- Alfvenic jets ^Ferreiral Il997|) . This has an important consequence 
on disc physics. Since i? + > B z , the total pressure gradient can overcome the 
magnetic compression due to B r only if /i is not larger than unity. The same 
constraint holds for the toroidal field which implies that a m must be of the 
order of unity. Finally, using the fact that A ^ r/h ^> 1 in a, JED, one obtains 
that fi cannot be too small. Thus, the parameter space for a JED is /i ~ q ~ 
a m ~ 1. Note also that a JED is thinner than a S AD because of the additio nal 
magnetic compression ijWardle fc KonigH Il993t IFerreira fc Pelletierl Il995|) . 
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Fig. 4. Magnetic acceleration arises whenever the projection of the Lorentz force on 
a poloidal field line becomes positive. This can be achieved in two ways, either with 
a downward vertical magnetic compression or a strong outward pressure force due 
to the toroidal field. The former leads to a small ejection efficiency and has current 
lines coming out of the disc surface (J z > 0) and entering at the inner radius. The 
latter has a strong ejection efficiency with the current entering the disc at it s surfaces 
(J z < 0). Only small ejection efficiencies allow for steady state solutions iFerreiral 



2.6 Disc radial equilibrium 

The quasi MHS radial equilibrium leads to an angular velocity 

dP/dr (J A B) r u 2 d\nu r \ 
pQ 2 K r pQ\r fl 2 K r 2 dlnr ) [ ' 

The deviation to the Keplerian rotation law Qk = \J GM/r 3 due to the radial 
(gas+radiation) pressure gradient is roughly of order (h/r) 2 at each altitude. 
This is because P scales with the density, which is not the case of the radial 
magnetic tension. At the disc midplane, it causes a deviation which is of the 
order ~ pR m (h/r) 2 ~ h/r but increasing as 1/p. Thus, thin accretion discs 
with h/r <C 1 will be mostly rotatin g at (sub-) Keplerian speeds but a problem 
arises when h ~ r (as in ADAFs iNaravan et alJlT998l or in self-gravitating 
discs, for instance). Indeed, it would imply a negative rhs at the disc surface 
which certainly means that no steady-state accretion-ejection solution can be 
found in that case. Note however that there is a priori no reason to ever have 
ft ~ r in JEDs: they are colder (see below) and more squeezed (see above) 
than a SAD. 



n 2 



n 2 



2.7 Energy budget 

The global energy budget is obtained by applying the energy conservation 
equation to the whole volume occupied by the JED. This equation writes 
P acc = 2P M hd + ZPrad where 

JiTi 
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is the mechanical power liberated by the accreting material between r e and 
r i, Pmhd — f Smhd ' dS is the flux through one disc surface of the MHD 
Poynting vector Smhd = E A B/fi a ~ — QxrB^Bp/ So, energy conser- 
vation in a JED tells us that the available accretion power is shared between 
a flux of electromagnetic energy powering the jets and radiation due to heat 
dissipation within the disc (2P ra d = Pdiss)- This dissipation is due to the fact 
that in a disc where turbulent magnetic diffusivity/resistivity and viscosity 
are assumed, there is always some heat production. The simplest and crud- 
est way to estimate this dissipation is using effective transport coefficients so 
that it writes Pdiss = J v dV^mJ'i+'q^Jp +r] v {rdfl / dr) 2 ), namely Joule and 
"viscous" heating. 

Since the "viscous" torque is negligible with respect to the jet torque, only 
a small fraction of the energy will be dissipated by viscosity. On the other 
hand, the most interesting aspect of angular momentum removal by jets is 
that the associated Joule dissipation implies also only a small fraction of the 
available ener gy. As a consequence, most of the liberated accretion power goes 
into the jets |Ferreira & Pelletieii fl993l Il995|) ! Precisely, this can be written 
as 

2P M hd A 2P rad 1 , 10 , 
~ and ~ - 12 

Pace 1 + A Pace 1+ A K ' 

where the ratio of the jet to the viscous torque A ~ r/h 3> 1. This property 
of JEDs has two important consequences: (i) the disc itself being weakly dis- 
sipative, it may well be unobservable leading to the (wrong) idea that there 
is no disc; (ii) a JED is cooler than a SAD fed with the same accretion rate, 
which leads to a smaller aspect ratio h/r. 



2.8 Links between jet and disc physics 

The previous sections showed the crucial role played by the magnetic diffusiv- 
ity within the turbulent JED. On the contrary, jets are best described by an 
ideal MHD formalism (y m = v' m = 0). This leads to the existence of 5 invari- 
ants along each magnetic surface for polytropic jets 2 . The Bernoulli equation 
is obtained by projecting the momentum equation @ along B p whereas the 
transfield or Grad-Shafranov equation by projecting it along Va (perpendic- 
ular to B p ). For more details see Tsinganos' contribution (this volume). 

MHD simulations of jets driven by accretion discs usually assume mag- 
netic field lines rotating at Keplerian speeds and negligible enthalpy leaving 
therefore 3 free and independent bou ndary conditions to be specified at each 
radius (see e.g. lAnderson et "aT]l2005l and references therein). These are often 
the density p(r), vertical velocity u z {r) and magnetic field B z (r) distributions. 
However, the study of MAES shows that not all distributions allow for steady 
state jets: there is a strong interplay between the disc and its jets. Such an 

2 The magnetic surface rotation rate fi(a), the mass to magnetic flux ratio ??(a), 
the total specific angular momentum L(a), energy E(a) and entropy K(a). 
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interplay appears in the form of analytical links between j et invariants and 
parame t ers descr i bing t he di sc. These links can be fo und in lCasse fc Ferreiral 
l|2000aft : iFerreiral 42002ft and lFerreira fc Cassd (|2004^ . Of all disc parameters 
the disc ejection efficiency £ plays a major role. Indeed, the knowledge of £ 
allows to define almost all jet properties. But in order to obtain the allowed 
values for £, the full set of MHD equations must be solved. 



3 A glimpse on self-similar solutions 
3.1 Mathematical method 

This is done by a separation method allowing to transform the set of partial 
differential equations (PDE) into two sets of ordinary differential equations 
(ODE) with singularities. Now, the gravitational potential in cylindrical co- 
ordinates is 

^ / x GM ( z 2 Y 1/2 

$a(r,z) = -— (^1 + ^J (13) 

and it is expected to be the leading energy source and force in accretion discs. 
Thus, if JEDs are settled on a large range of radii (so that we do not care about 
the radial inner and outer boundaries) , then the magnetic energy density has 
to follow gravity in order to match it everywhere. It is therefore justified to 
look for solutions of the form A(r, z) = GA{r)fA(f:) for any physical quantity 
A(r, z). Moreover, since gravity is a power law of the disc radius, we will use 

the self-similar Ansatz A(r,z) = A e I — j Ia{x) where x = zjh(r) is our 
self-similar variable with h oc r and r e is the JED outer radius. Because all 
quantities have power law dependencies, the resolution of the "radial" set of 
equations is trivial and provides algebraic relations between all exponents. 
The most general set of radial exponents allowing to take into account all 
terms in the dynamical equations leads to the following important constraint 

0=1 + 1 (14) 

where the magnetic flux distribution writes a(r) o c r@ . As an illustration, 
the solutions obtained bv lBlandford fc Pavnel 1 1982ft used (3 = 3/4, ie £ = 0. 
In general, all self-similar models of disc driven jets not addressing the disc 
dyna mics use a magnetic field distribution inconsistent with the jet mass load- 



ing ; (Blandford fc Pavn 2 Il982t IContopoulos fc Lovelacel Il993 IContopoulosl 
ll994UOstrikeiill997tlVlahakis et all 1200(1) . 

All quantities Ja{x) are obtained by solving a system of ODE which can 
be put into the form 



d£i 

dx 

M 

dx 
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where M is a 8x8 matrix in resistive MHD regime, 6x6 in ideal MHD 
l|Ferreira k Pelletierl Il995|) . A solution is therefore available whenever the 
matrix M is inversible, namely its determinant is non-zero. Starting in resis- 
tive MHD regime, det M = whenever V 2 (V 2 — Cg) — where C s is the sound 

2 

speed and V = u ■ n is the critical velocity. The vector n — (e z — -e r )(l + ^) 
provides the direction of propagation of the only waves consistent with an 
axisymmetric, self-similar description (see Tsinganos' contribution). There- 
fore, close to the disc, the critical velocity is V ~ u z , whereas far from 
the disc it becomes V ~ u r (no critical point in the azimuthal direction). 
Inside the resistive disc, the anonalous magnetic resistivity produces such 
a dissipation that the magnetic force does not act as a restoring force and 
the only relevant waves are sonic. Note also that the equatorial plane where 

V = is also a critical point (of nodal type since all the solutions must pass 
through it). This introduces a small difficulty as one must start the integra- 
tion slightly above z = 0. In the ideal MHD region, defh/l — whenever 
(V 2 ~ v smW 2 ~ v fm)( v2 ~ V Anf = namely, when the flow velocity 

V successively reaches the three phase speeds Vsm, Vau and Vfm, corre- 
sponding respectively to the slow magnetosonic (SM), Alfven and fast mag- 
netosonic (FM) waves. The phase speeds of the two magnetosonic modes are 

V§m,fm - I (C 2 + VI T V(C 2 + Vl t f - 4CJV£) where V At is the total 
Alfven speed and Vau = V Ap ■ n. These expressions are slightly modified by 
the self-similar ansatz. Note however that the condition V — Vau is equivalent 
to u p = Va p - 




Fig. 5. Once material has left the resistive MHD zone, it is frozen in a particular field 
line and encounters the three MHD critical points. The smooth transition between 
resistive and ideal MHD regimes already selects the MAES parameter space (see 
Section 2). 



How do we proceed ? We fix the values of the four disc parameters 
(e = h/r,a m = v m /V A h,Xm = v m jv' m ,V m = v v /v m ) and some guesses for 
the disc magnetization fi and ejection efficiency £. Starting slightly above the 
disc midplane where all quantities are now known, we propagate the resistive 
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set of equations using a Stoer-Burlisch solver for stiff equations. As x = z/h 
increases, the flow reaches an ideal MHD regime and we shift to the corre- 
sponding set of equations. Care must be taken in order to not introduce jumps 
in the solution while doing so. The smooth crossing of the SM point can only 
be done with a critical value for /x. We thus modify our initial guess until the 
solution gets close enough to the critical point and jump across it (leapfrog 
method). The same must be done for the Alfven point which requires a criti- 
cal value for £. Each time another guess for £ is made, one has to find again 
the corresponding critical value for (x (Fig. 0) . The crossing of the last critical 
point (FM) does not bring much more information on MAES physics and will 
be discussed in Section 4. 

3.2 Typical solutions 

Only the most salient features of sel f-si milar accretion-ejection solutions will 
be discussed here (see lFerreirall2002l and iFerreira fe Cassdl2004|) . 

Cold solutions 

Cold solutions are defined h ere by an isothermal IFerreira fc Pelletierl Il995t 
lFerreiralll997|) or adiabatic ijCasse fc Ferreira! l2000a|) energy equation Since 
the plasma pressure P is (h/r) 2 smaller than the gravitational energy density, 
such an energy equation ensure s that the jet enthalpy is ne gligible with respect 
to gravity and magnetic fields llBlandford fc Pavnelllflsi . 

Figure [S] shows the velocity components in both the JED and the jets as 
a function of the self-similar variable x, along a magnetic surface for typical 
solutions with h/r — 0.01 but different ejection efficiencies £. The disc surface 
is located at x = 1 and the Alfven point is reached at x ~ 100 (za ~ r/Cj. 
Note that the disc vertical velocity is negative within the disc (material is 
falling) and becomes positive only slightly before the point where the radial 
velocity itself becomes positive. This happens roughly at the disc surface (see 
bottom right panel) but still in the resistive MHD regime. The SM point is 
crossed at x ~ 1.6. All velocity components are comparable at the Alfven point 
(this also holds for the magnetic field). Beyond that point, the plasma inertia 
overcomes the magnetic tension and the magnetic surface opens tremendously. 
This leads to the build up of a sheared magnetic configuration (the ratio 
\B<j,/Bp\ increases). Note that the structure of the jet can be characterized 
by two families of intertwined helices: the plasma streamlines (wound in the 
same direction as the disc rotation) and the magnetic field lines (wound in 
the opposite direction). 

The magnetic acceleration is so efficient that all available MHD energy 
is transferred into jet kinetic energy. From Bernoulli equation one gets an- 
alytically the asymptotic jet velocity Vj — f2 r y/2A — 3 where f2 r is the 
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Fig. 6. Components of the jet poloidal velocity u v and logarithm of the ratio 
of the poloidal to the azimuthal velocity, measured along a magnetic surface for 
£ = 0.005 (solid line), 0.01 (dotted line), 0.02 (short-dashed line) and 0.05 (long- 
dashed line) (e = 10 -2 , a m = 1). For these typical cold solutions, the jet always 
reaches its maximum velocity, mainly as a vertical component (the jet opening angle 
is tan# = B r /Bz = u r /u z ). Inside the disc, matter is b eing accreted with a velocity 
of order e the Keplerian velocity Q r l|Ferreiral 11997"). 

Kepleria n speed at the jet footpoin t r a and A is the magnetic lever arm pa- 
rameter (jBlandford fc Pavnel Il982() . This importan t jet parameter is actually 
related to the disc ejection efficiency A ~ 1 + l/2£ l|Ferreiralll997jl . 

The disc parameter space has been thoroughly investigated for cold solu- 
tions. It is very narrow with typical values £ ~ 0.01 and 0.1 < fj, < 1, with 
the following approximate scaling 

£~0.1 M 3 (15) 

Although its validity holds only in a quite narrow interval, it shows that the 
stronger the field the more mass is ejected. No solution has been found outside 
the range 0.0007 < e = h/r < 0.3 and 0.3 < a m < 3 (FigEJ. As pointed out 
previously, there is no solution with a dominant viscous torque. All solutions 
exhibit a high degree of collimation: actually, they even undergo re collimation 
towards the axis which should result in a shock l|Ferreiral Il997|) . However, 
the subsequent behaviour of the jet after that shock cannot be treated within 
self-similarity. 
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Fig. 7. Disc parameter space for is othermal jets j Ferreir al Il997l) (adiabatic jets 
present no qualitative difference, see lCasse fc Ferrejrell^nnnah . Left: disc magneti- 
zation (i as a function of the disc ejection efficiency £ for a m — 1 and various disc 
aspect ratios: e = h/r — 10 _1 (solid line), 10 -2 (dotted line), 10 -3 (short-dashed 
line) and 7 10 -4 (long-dashed line). The main effect of decreasing e is to shift the 
range of allowed £ to higher values (but with a more limited range). Right: Influence 
of the turbulence parameter a m on the disc magnetization fi for e = 1CP 1 and vari- 
ous ejection efficiencies: £ = 0.004 (solid line), 0.005 (dotted line), 0.01 (short-dashed 
line) and 0.02 (long-dashed line). The minimum level of MHD turbulence is limited 
by the value of the induced toroidal field allowing trans- Alfvenic jets, whereas the 
maximum level has been arbitrarily fixed to unity. 



Warm solutions 

Warm solutions are obtained by solving Eq. 10 with a prescribed self-similar 
function Q. Several physical effects can be simulated that way: 

• Heat deposition at the disc surface only: the function Q reaches a maximum 
at the disc upper layers and then decreases rapidly (to recover adiabatic 
jets). This mimics the effect of disc illumination by stellar UV and X 
rays. Alternatively, this energy could arise from the dissipation of a small 
fraction of the accretion energy, released in these layers by turbulence. 
Remarkably the mass loa d can be significantly en hanced, with ejection 
efficiencies up to £ ~ 0.46 ()Casse fc FerreiraLl2000rl . 

• Heating of the sub-Alfvenic regions: the function Q is non zero in these 
regions only with subsequent adiabatic or polytropic jets. This mimics 
the effect of some "coronal" heating as in the solar wind or, alternatively, 
the pressure due to an inner flow (e.g. stellar or magnetospheric wind) 
ramming into the disc wind. Under some circumstances, the field lines are 
forced to open much more than they would which results in a different jet 
dynamical behaviour. In particul ar, self-similar je t s can smoothly cross the 
last m odified FM critical point (|Vlahakis et all l200ft iFerreira fc CasseL 
l2004h . See Fig. ED for an example. 

• Heating of the whole jet: this has not yet been done in the framework of 
disc driven jets (but could easily be done by assuming a positive function 
Q everywhere in the jet). The reason is that such jets would not be sig- 
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Fig. 8. Typical super-FM disc wind with £ = 0.03, e = 0.03 (ft = er). Density, 
pressure and temperature are normalized to their value at the disc midplane, the 
magnetic field components to B z (z = 0) and the velocities to the Keplerian speed 
at the anchoring radius r a . All magnetic field components remain comparable from 
the disk surface to the Alfven point. Note that the density profile inside the disc, 
where both u T and u z are negative, is very different from a gaussian. Recollimation 
takes place at z ~ 3 10 3 r o . The lower right panel shows the various critical Mach 
numbers (e.g. Msm = V/Vsm) appearing in the self-similar equations. The usual 
fast Mach number, n = u v /u fm, becomes greater th an unity much sooner than the 
critical one Mfm = V/Vfm l|Ferreira fc CasselEooih . 



nificantly modified by a warmer material (in contrast with stellar winds) . 
However, this is interesting for comparing models to observations. Indeed, 
observed jets display temperatures of some 10 4 K that require some heat- 
ing mechanism(s) overcoming the huge cooling due to the jet expansion 
(so called adiabatic cooling). It has been shown that ambipolar diffusion 
is not enough and t hat some turbulent or shock heating must be at work 
l|Garcia et alll200l[) . 



4 Concluding remarks 
4.1 What's next ? 

The theory of steady jet production from Keplerian accretion discs is now com- 
pleted. The physical conditions required to thermo-magnetically drive jets are 
known, all relevant physical processes have been included in the framework of 
mean field dynamics. Of course, there are still many unsolved questions: 
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(i) Can a sustained MHD turbulence maintain a m ~ 1? This is a huge con- 
straint that deserves a thorough investigation. 

(ii) Obse rvations of T Tauri jets favor solutions with large ejection efficiencies 
(£ ~ O.l. lPesenti et al.ll2004j) requiring additional heating at the disc surface. 
A theoretical assessment of this heating must be undertaken. 

(Hi) What is the stability of MAES? As will be seen below, there was some 
claims that MAES were unstable but they were proven to be wrong. On the 
other hand, jets do show time dependent features and one must clearly go 
beyond steady state models. On that respect, numerical simulations will be 
very helpful. 

(iv) Disc driven winds do not treat the star-disc interaction. Understanding 
the whole process of star formation requires now to address this crucial issue 
as it pinpoints the problem of the stellar angular momentum removal. This is 
further discussed below. 



4.2 Biases of self-similarity 

Self-similarity allows to take into account all dynamical terms in the equations 
and, as such, is the best means to solve in a self-consistent way the steady-state 
accretion-ejection problem. However, there is a price to pay... 

(i) The asymptotic behaviour is obviously biased since, for instance, 
neither inner nor outer pressures can be taken into account. In fact, no realistic 
"radial" boundary condition whatsoever can be dealt with. When modeling 
an astrophysical jet, this implies for instance to truncate the solution at one 
inner and outer radius. But there is anoth er aspect , less known and more 
subtle. IContopoulos fc Lovelacel obtained jet so- 
lutions within the same self-similar framework but with different asymptotic 
behaviors. The reason stems from the fact that they played around with j3 
(flux function a(r) oc r") as if it were a free function whereas the mathemat- 
ical matching with a Kep lerian disc imposes its value. On the other hand, 
IPelletier fc Pudrit3 l|l992|) obtained also recollimating non self-similar solu- 
tions, which indicates that recollimation can indeed be physical and not en- 
tirely due to self-similarity. In fact, it can be shown that recollimation of a 
jet launched from a Keplerian accretion disc is poss ible whenever the radial 
profile of the ejection efficiency £ is smooth enough <)FerreiraLll997l) . 

(ii) The regularity conditions are to be imposed at the modified points 
and not at the usual ones (see Tsinganos, this volume). However, these loca- 
tions coincide for both the slow (SM) and the Alfven points so that one can 
be confident that there is no bias there. However, this is not so for the fast 
magnetosonic point. Se lf-similar trans-FM solut ions require an Alfven surface 
very close to the disc ()Ferreira fc CasseL 120041) , which can only be done by 
the action of a large pressure in the sub-Alfvenic region. This is obviously a 
strong bias since it is not clear whether such a pressure is indeed provided 
in astrophysical objects. Note however that crossing this modified FM point 
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is more a theoretician satisfaction than anything else: it gives no additional 
physical insight on the disc physics. 

(iii) The local disc physical conditions as obtained with self-similar 
solutions are not biased. The physical processes are well identified and un- 
derstood and can be sometimes even obtained in a pure analytical manner. 
They have b e en also confirmed b y numerical experiments of lCasse fc Kennensl 
l|2002l 12004ft : IZanni et all l|2004l) (although one might object that numerical 
experiments were actually tested with the help of semi-analytical solutions). 

4.3 Is accretion-ejection unstable? 

There has been some cl aims in the literature that the accretion-ej ection pro- 
cess itself was unstable l|Lubow et all Il994bl ICao fc Spruil 1200211 . The idea 
was the following. Start from a steady picture where the accretion velocity 
u r at the disc midplane is due to the jet torque. It leads to a bending of the 
poloidal field lines described by an angle 9 with the vertical. Now imagine a 
small perturbation Su r enhancing the accretion velocity. Then, according to 
these authors, the field lines would be more bent (9 increases) which would 
lead to lower the altitude of the sonic point. Because the sonic point would 
be located deeper in the disc atmosphere, where the density is higher, more 
mass would be henceforth ejected which would then increase the total angular 
momentum carried away by the jet. This means that the torque due to the 
jet is enhanced and will, in turn, act to increase the accretion velocity. Thus, 
the accretion-ejection process is inherently unstable. 

The whole idea of this instability is based upon a crude approximation 
of the disc vertical equilibrium. In fact, the magnetic field produces a strong 
vertical compression so that, as 9 is i ncreased, less mas s is be ing ejected, 
not m ore. This has been pointed out bv lKonigl fc Wardlel l)l996j) and iKonigl 
l)2004|) and is indeed verified in the full MAES calculations reported here. 

4.4 Magnetic fields in accretion discs 

The necessary condition for launching a self-collimated jet from a Keplerian 
accretion disc is the presence of a large scale vertical magnetic field close to 
cquipartition l|Ferreira fc Pelletierl [l995) , namely 

This value is far smaller than the one estimated from the interstellar magnetic 
field (see Menard's contribution), assum i ng either id eal MHD or B oc n 1 / 2 
i Heiles et all Il993t ISasu fc MouschoviasU 19941 119951) . This implies some de- 
coupling between the infalling/ accreting material and the magnetic field in 
order to get rid off this field. This issue is still under debate. The question is 
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(a) Extended disc-wind: r e »rj (b) X-wind: ie>i"i 




Fig. 9. Two classes of stationary accretion powered disc winds, (a) "extended disc 
winds", when the magnetic flux threading the disc is large enough so that a large 
radial extension of the whole accretion disc drives jets (r e S> r*i). The Alfven surface 
is expected to adopt a rather conical shape, (b) "X- winds", when the magnetic 
flux is small and only a tiny disc region is driving jets. The Alfven surface can be 
either convex or concave, although the latter is probably more physical (since less 
material can be ejecte d at the two extremes and the Alfven point is rejected to 
infinity). Adapted from iFerreira et alJ (I2006al) . 

therefore whether accretion discs can build up their own large scale magnetic 
field (dynamo) or if they can drag in and amplify the interstellar magnetic 
field? Although no large scale fields have been provided by a self-consistent 
disc dynamo, this scenario cannot be excluded. But the latter scenario (ad- 
vection) seems a bit more natural. 

Let us assume that the disc material is always ionized enough to allow for 
some coupling with the magnetic field (and use MHD). The outer parts of the 
accretion disc will probably t ake the form of a SA D with no jets and almost 
straight (TZ m ~ 1) field lines l|Lubow et all Il994al) . In that case, the steady- 
state solution of the in duction equation for the poloidal field is B z oc r~' Rm 
ijFerreira et all l2006bj) . Hence, as a result of both advection and (turbulent) 
diffusion, the magnetic field in a SAD will be a power-law of the radius. 

Can a SAD transport B z and allow for a transition to an inner JED? 
This will be so if there is some transition radius (the outer JED radius r e ) 
where /i = B 2 z j /i D P becomes of order unity In a SAD the total pressure writes 
p _ MaOjch K r -3/2-(5 w j^ n ^ r s _ g mce g j s a l wa y S close to unity in 

circumstellar discs, one gets /i a r f with e ~ 1. Thus, it can be readily 
seen that it is indeed reasonable to expect such a transition (computing it 
is another matter), at least in some objects. The recent Zeeman observation 
of a magnetic field in the accretion disc of FU Or supports this conclusion 
l|Donati et all 12005^1 . 

4.5 X-winds and disc winds 

The X-wind model llShu et all Il994t iNaiita fc Shil Il994l IShu et all 1199.1: 
IShang et al ., 1998. 120021) is a rich and complex model but, contrary to common 
belief, it is an accretion-powered wind launched from the accretion disc. In 
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practice, if the amount of magnetic flux threading the disc is large so that 
r e >• then one gets an "extended disc wind", whereas if the magnetic 
flux is tiny with r e > r i: one gets an "X-wind" (Fig. |5J). The dynamics and 
asymptotic behaviour of jets will differ strongly between an ext ended disc wind 
and an X-wind and can thereby be tested against observations l|Ferreira et all 
2006a). But this difference arises mainly because of the restricted range in 
radii in the X-wind case, not because the underlying disc physics is different. 
The basic phenomena described in Section 2 apply as well for the portion 
of the disc launching the X-wind. Thus, equipartition fields are required, the 
"viscous" torque is negligible with respect to the jet torque and the angular 
momentum carried away by the X-wind is exactly the same amount lost by the 
accreting material. As a consequence, X-winds cannot take away any angular 
momentum from the central star. 

Published material on the dynamics of X-winds contains: (i) a scenario for 
the origin of B z (stellar) and the star-disc interaction (leading to the opening 
of some magnetic field lines); (ii) the calculation of the sub-Alfvenic ideal 
MHD jet (elliptic domain defined by prescribed boundary conditions); (iii) a 
somewhat mysterious "interpolation" to a simple jet asymptotic solution. The 
following questions remain therefore to be addressed: 

(1) Can the disc afford the imposed mass flux and field geometry? Indeed, 
the assumed ejection to accretion mass flux ratio of 1/3 from such a tiny 
region is huge and would require a fantastic ejection efficiency (£ of order 
unity or larger). The calculations of JEDs showed that this is unfeasible in 
a steady way. However, the huge magnetic field gradients required in the 
X-wind launching region provide a significantly different situation. This has 
never been analyzed. 

(2) How good is the transfield equilibrium satisfied? There is no mathematical 
procedure to find a solution of mixed type (elliptic-hyperbolic) PDEs when 
the singular surfaces are unknown. The trick used for X-winds provides an 
incomplete solution, but there is maybe some means to fulfill the transverse 
equilibrium by using an iterative scheme. In any case, this important point is 
missing in the current published material. 

4.6 Magnetic star-disc interactions 

Nowadays it seems accepted that a lot if not all young stars have a mag- 
netospheric interaction with their circumstellar accretion disc (see Alencar's 
contribution, this volume) . If one assumes that the disc is threaded by a large 
scale magnetic field, then the question of how this field is connected to the 
stellar field arises. First ideas are always simple and so is the stellar magnetic 
field, assumed up to now to be dipolar and axisymmetric (see Mohanty et 
al. 2006). We define here the magnetopause as the radius r m below which all 
field lines threading the disc are tied to the star whereas beyond r m , they are 
disconnected from the star. 
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(a) 



Y-type Interaction 



(b) 



X-type Interaction 




Fig. 10. Two simple axisymmetric star-disc magnetospheric interactions, (a) "Y- 
type" interaction obtained when the stellar magnetic moment is anti-parallel to the 
disc magnetic field. A current sheet is formed at the interface between the open 
stellar field and the disc field. Such a configuration cannot produce per se a wind, 
(b) "X-type" configuration obtained when the stellar magnetic moment is parallel 
to the disc magnetic field. A magnetic X-point is generated at the disc midplane 
where the two fields cancel each other. Unsteady e jection (" ReX- winds" ) can be 



The case envisioned within the X-wind scenar io assumes a stell ar magnetic 
moment anti-parallel to the disc magnetic field l|Shu et all Il994j) . As a con- 
sequence, a neutral surface (where B = 0) appears above each disc surface, 
illustrated by a limiting poloidal field with a Y shape (Fig. llOfl . The other 
case, a stella r magnetic moment p arallel to the disc magnetic field, has been 
proposed bv iFerreira et all (|200Cf) . The two fields then cancel each other at 
the disc midplane, defining a neutral line at a radius rx where reconnection 
takes place. This configuration gives rise to "Reconnection X- winds" (here- 
after ReX-winds) specifically above this zone. 

Accretion curtains 

The first question here is can these simple topologies allow for accretion be- 
low 3 r m ? Disc material will accrete only if it looses angular momentum and 
this depends on both turbulence and the magnetic torque due to the mag- 
netosphere. The magnetosphere will try to make the disc corotate with the 
protostar so the sign of the torque depends on their relative angular veloc- 
ity. The corotation radius, r co — (GM/^) 1 '' 3 , is defined as the radius where 
the stellar angular velocity is equal to the Keplerian one. This gives an 
estimate of the real angular velocity of the disc (since the disc magnetic field 
introduces already a deviation) . Roughly speaking, if r m > r co the star rotates 
faster than the disc and deposits its angular momentum, whereas if r m < r co , 
the star rotates slower and thus spins down the disc. Note that r m denotes 
roughly the radius where the stellar magnetic field becomes dynamically dom- 



3 Accretion is realized beyond r m by e.g. the jet torque within the JED and, farther 
away in the SAD, by the turbulent "viscous" torque. 
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inant, namely /z > 1. Thus, unless a very efficient turbulent mechanism is 
operating and transports radially the stellar angular momentum, no accretion 
is possible when r m > r co (although such a "propeller" regime is favorable for 
ejection). 

As a consequence, both X-type and Y-type interactions allow for a mag- 
netospheric accretion as long as r m < r co . It is interesting to note that both 
configurations require an equatoria l reconnection zone (interesting for sud den 
energy dissipation and chondrules, IShu et aljEnmllnourielle et al.ll'.>()()(ili. In 
the case of a Y-type interaction, it arises because of the requirement that the 
magnetospheric field makes an angle with the vertical large enough in order to 
allow the disc material to flow inwards. This assu mption implies a mag netic 
neutral "belt" at the disc midplane (see Fig. 1 in lOstriker fc IiE]|l995), but 
whose origin and dynamics were not discussed and remain therefore major 
unsolved issues. In the case of an X-type interaction, the presence of the mag- 
netic neutral line is due to and maintained by the cancellation of the two fields 
(see fig llOf) . The accreting disc material can cross the resistive MHD region 
and is lifted vertically by the strong Lorentz force above the reconnection 
site. The transition from an accretion disc to accretion curtains can be quite 
smooth in that case. 

Stellar spin down 

The second q uestion is the iss u e of t he stellar angular momentum removal 
by winds (see iMatt fc Pudritd l|2005j) for more details and the necessity of 
winds). As explained earlier, X- winds carry away the angular momentum of 
the accreting disc material. Thus, such a configuration cannot brake down 
the protostar (as initially claimed). On the contrary, the X- type configuration 
provides a very efficient means to do it l|Ferreira et allEoOOj) . The reason is the 
possibility to launch disc material above the reconnection site. The scenario is 
the following (Fig. lll|) . A stationary extended JED is settled in the innermost 
regions of the accretion disc and provides open magnetic flux to the star. This 
magnetic field reconnects at rx with closed stellar field lines: the disc field 
contributes thereby to transform closed magnetospheric flux into open flux. 
At the reconnection site, the disc material is lifted vertically and loaded onto 
these newly opened field lines, tied to the rotating star. Whenever rx > r co , 
the star is rotating faster than the loaded material and it undergoes a strong 
magneto-centrifugal acceleration. This gives rise to the so-called ReX-wind, 
whose energy and angular mo mentum are th o se of t he star. Using a toy-model 
for the magnetic interaction iFerreira et alJ l)2000(l showed that such winds 
could brake down a contracting protostar on time scales that are comparable 
to the duration of the embedded phase (Class and I sources). The protostar 
was assumed to rotate initially at breakup speed and, after some 10 5 to 10 6 

4 Note that it should be operating whe n fx > 1, while the magneto-rotational in- 
stability is already quenched at [i ~ 1 fealbusU2003f) . 
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Cylindrical radius R (AU) 

Fig. 11. The ReX-wind configuration jFerreira et alll20odl . A MAES is established 
around a protostar whose magnetic moment is parallel to the disc magnetic field. 
This is a natural situation if both fields (disc and stellar) have the same origin. 
Left: black solid lines are streamlines, white dashed lines are contours of equal total 
velocity (mainly rotation inside the disc) and the background color scale shows the 
density stratification. The ReX-wind (arrows) would be confined and channeled by 
the outer disc wind. Right: sketch of the magnetic configuration leading to Rex- 
winds and accretion curtains around the magnetic neutral line at rx- Arrows show 
the expected time-dependent plasma motion. 

yrs, it has been spun down to 10% of it despite its contraction and mass 
accretion. 

ReX-winds seem therefore to offer a serious possibility to brake down pro- 
tostars (to my knowledge, there is no other model in the literature). Note that 
ReX-winds are probably intermittent by nature because of the unavoidable 
radial drift of the reconnection site (there is no ejection whenever rx < f co ). 
Dynamically speaking, such an unsteady "wind" should be better described 
as bullets flowing inside the hollow disc wind. Remarkably the basic features 
of X-type configurations remain if the stellar dipole is inclined: one would ob- 
serve in that case precessing bullets channeled by the outer disc wind. Heavy 
numerical simulations will be required to test and analyze this scenario. 
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